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摘要 受 地 球 大 气 影响 , 使 用 地 基 光 学 望远镜 观测 角 距离 较 小 的 双星 系统 或 主 卫星 系统 时 往往 会 出 现 星象 不 
可 分 辩 的 情况 . 因此 , 系统 光 心 位 置 与 系统 质心 位 置 可 能 存在 一 定 的 偏 莽 . 准确 地 测量 太阳 系 天 体系 统 质心 位 
置 可 以 改进 其 轨道 参数 , 有 助 于 揭示 太阳 系 的 形成 与 演化 . 以 矮 行 星 Haumea 及 其 亮 卫 星 Hi'iaka 的 运动 为 例 , 仿 
真 系统 光 心 围绕 质心 扰动 的 过 程 , 探 完 视 宁 度 ( 用 星象 的 半 高 全 宽 表 示 ) 变 化 对 准确 测量 光 心 位 置 的 影响 . 仿真 
结果 表明 , 使 用 二 维 高 斯 定 心 算法 测定 的 系统 光 心 位 置 随 视 宁 度 变化 , 而 修正 矩 定 心 算法 的 定 心 结 果 不 受 视 
宁 度 的 影响 . 根据 仿真 结果 , 研究 能 够 有 效 减 少 视 宁 度 变 化 对 交心 位 置 准 确 测量 影响 的 定 心 算 法 十 分 必要 ; 同 
时 , 新 的 定 心 算法 还 需 考虑 主星 光度 变化 的 影响 . 使 用 云南 天 文 台 2.4 m 望 远 镜 , 1 m 望 远 镜 以 及 紫金 山 天 文 台 
姚 安 观测 站 0.8 m 望 远 镜 从 2022 年 2 月 7 日 至 2022 年 5 月 25 日 观测 矮 行 星 Haumea 系 统 , 得 到 29 晚 共 463 幅 CCD 图 
像 . 新 定 心 算 法 确定 的 光 心 位 置 与 使 用 二 维 高 斯 定 心 算法 的 结果 相 比 具有 更 好 的 位 置 拟 合 效 末 . 此 外 , 还 发 
现 亮 卫星 Hi'iaka 在 Jet Propulsion Laboratory (JPL) 历 表 与 Institut de Mécanique Céleste et de Calcul des 


Éphémérides (IMCCE) 历 表 中 的 位 置 存在 较 大 偏差 . 


关键 词 天 体 测量 学 , 技术 : 图 像 处 理 , 方法 : 数据 分 析 , 柯 伊 伯 带 天 体 : 个 别 : 妊 神 星 


中 图 分 类 号 : P129;， 文 献 标识 码 : A 


1 引言 

位 于 太阳 系 海王 星 轨道 之 外 的 柯 伊 伯 币 结构 
巨大 , 其 研究 对 于 揭示 太阳 系 的 形成 和 演化 具有 
重要 意义 . 许多 柯 伊 伯 带 天 体 拥 有 伴星 或 卫星 (以 
下 简称 主 卫星 系统 ), 但 当 使 用 地 基 光 学 望远镜 观 
测 它们 时 , 其 在 CCD 图 像 上 往往 表现 为 单一 的 混 
合 星象 , 比如 经 过 国际 天 文学 联合 会 (International 
Astronomical Union, IAU) WE HI & E Æ (Pluto), 
4H Æ (Makemake), 14 Œ (Haumea) 7847 42 
其 卫星 构成 的 主 卫 星系 统 都 存在 此 种 现象 . 即 在 
A} HF EAS EIN, 观测 得 到 的 是 系统 光 心 位 置 . 当 

2023-08-28 收 到 原稿 , 2023-12-11 收 到 修改 稿 


主星 和 卫星 的 质量 比 与 光度 比 相 差 较 大 时 , 系统 
光 心 和 系统 质心 将 出 现 不 重合 的 情况 . 197048, 
Morel 岂 表明 光 心 绕 转 质心 作 开 普 勒 椭圆 运动 . 
2011 年 ，Ortiz 等 站 在 对 大 柯 伊 伯 带 天 体 Orcus (有 
一 颗 卫 星 Vantb) 进 行 33 qd 的 地 面 观测 中 发 现 了 其 系 
统 光 心 位 置 存 在 周期 扰动 现象 , 即 晃 动 (wobble). 
通过 长 期 高 精度 天 体 测量 得 到 的 准确 光 心 运动 轨 
迹 在 双星 系统 中 常 被 用 于 双星 轨道 求解 、 质 量 计 
算 中 以 及 实际 质心 位 置 推 导向 , 该 方法 同样 适用 于 
主 卫 星系 统 . 

准确 地 测量 系统 光 心 位 置 依赖 于 许多 因素 . 据 
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van de Kamp, 系统 光 心 在 理论 上 是 主星 和 
卫星 光度 的 加 权 中 心 点 , 光 心 到 主星 的 角 距 离 xr1 
只 取决 于 主 卫 星之 间 的 星 等 差 及 角 距 离 . 但 从 地 
面 观测 时 , 光 心 的 准确 位 置 测 量 可 能 受到 大 气 视 
宁 度 变化 的 影响 . 本 文通 过 计算 机 仿真 方法 模拟 
了 主 卫 星系 统 运动 过 程 , 详细 探究 系统 光 心 位 置 
测量 受 视 宁 度 变 化 的 影响 及 相应 的 校正 方法 . 本 
文选 用 矮 行 星 Haumea 及 其 完 卫 星 Hi"iaka 构 成 的 主 
卫星 系统 作为 研究 对 象 . 由 于 Haumea 另 一 箱 暗 卫 
星 Namaka 的 视 星 等 约 为 21.8 magl, 超过 了 本 次 
探究 所 使 用 望远镜 CCD 在 适度 曝光 时 间 内 的 极限 
ORR S), 本 文 忽 略 了 其 带 来 的 位 置 测量 影响 . 
矮 行 星 Haumea 系 统 是 柯 伊 伯 带 合 撞 家 族 中 已 知 的 
mA BYP), 所 以 主星 Haumea 和 卫星 Hi'iaka 应 具 
有 相似 的 物理 性 质 , 它们 表面 宪 盖 着 大 面积 的 水 


Je RUU 主星 Haumea 及 其 卫星 Hi'iaka 的 平均 星 
等 差 为 2.76 mag (J、H 波 段 ). Haumea 系 统 的 反照 
率 数值 接近 于 1 且 为 碰撞 家 族 , 表明 了 主星 和 卫 
星 颜 色相 近 ( 在 可 见 光 和 近 红 外 波段 ). 因此 , 它们 
之 间 的 星 等 差 是 相近 的 . 本 文采 用 了 平均 星 等 差 
为 2.76 mag ( 即 主 卫星 通 量 比 为 12.71) 进 行 试验 . 卫 
星 Hi"iaka 绕 主星 的 公转 周期 为 (49.12 + 0.03) d}, 
两 者 之 间 的 角 距 离 在 0.44” -1.45” 之 间 变 化 ( 取 自 Jet 
Propulsion Laboratory (JPL) 历 表 ! 和 Institut de 
Mécanique Céleste et de Calcul des Ephémérides 
(IMCCE) 4"). 主星 与 亮 卫星 Hi'iaka 的 质量 比 为 
223.80 + 13.7504. 此 外 , 矮 行 星 Haumea 有 明显 双 
峰 自转 光 变 曲线 9, 主星 自转 引起 的 光度 变化 对 系 
统 光 心 的 影响 也 应 该 考虑 在 内 . 

本 文 其 余部 分 组 织 结构 如 下 . 第 2 市 将 阐述 主 
卫星 系统 运动 的 规律 , 基于 仿真 实验 的 结束 提出 
校正 光 心 的 新 定 心 算法 ; 第 3 节 将 介绍 Haumea 系 
统 的 观测 资料 和 归 算 过 程 , 得 出 光 心 校正 后 的 位 
置 测量 结果 ; FATS on Hr AT VO BS IE BR K E 
卫星 Hi'iaka 在 JPEL 历 表 和 IMCCE 历 表 的 位 置 差异 ; 
最 后 一 节 进 行 总 结 与 展望 . 


"https: //ssd.jpl.nasa.gov/horizons.cgi#top / 


*http://vo.imcce.fr /webservices/miriade/ 
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2 ”研究 方法 
2.1 ” 主 卫星 系统 概述 

图 1 描绘 了 主 卫 星系 统 的 运动 过 程 . 假设 系 
统 的 质心 Cu 固定 不 动 , 主星 P (Primary) 和 卫星 S 
(Satellite) 同 时 绕 质 心 做 同一 方 同 的 开 普 勒 椭圆 运 
动 , 光 心 Cr (五 角 星 形 ) 绕 质心 Cv ( 圆 形 ) 也 做 相同 
方 回 的 椭圆 运动 . 图 中 , 光 心 和 质心 到 主星 的 角 距 
离 分 别 表 示 为 已 和 rw 光 心 到 质心 的 角 距 离 表示 为 
A, 主 卫 星之 间 的 角 距 离 表 示 为 a. 据 文 献 [5], 光 心 
位 置 由 主星 和 卫星 的 光度 共同 决定 , 始终 在 两 者 
的 连 线 上 , 更 靠近 光度 较 大 的 主星 (符合 杠杆 原理 ). 
因此 如 (1) 式 ， rp ake th, 其 中 Lp 和 Ls 分 别 为 主 
星 和 卫星 的 光度 . 


1 
=e, (1) 
1+ 


rv 遵循 相同 原理 , 则 A 的 表达 式 如 下 , 其 中 Mp 和 
Ms 分 别 为 主星 和 卫星 的 质量 
a | 
(+75) (+456) 


TL 


= ha. 


(2) 
ECHR], oe - hp RAMA. 当 
ABA EIT, et SOS Pk) HOH 
圆 ) 形 状 与 卫星 绕 主 星 公转 轨迹 (最 大 的 椭圆 ) 形 状 
相同 , 两 者 的 大 小 成 比例 . 从 地 面 观测 , 唯 有 系统 光 
心 CL 是 视觉 上 可 见 且 实 践 中 可 测量 的 , 所 以 图 1 也 
解释 了 在 对 某 些 目标 长 时 间 观 测 时 , 观测 的 光 心 位 
置 会 出 现 周期 性 波动 现象 . 


2.2 ” 光 心 与 视 宁 度 变 化 关系 探究 

本 文 沿用 Ortiz 等 习 的 仿真 方法 , 使 用 photutils 
库 3 模 拟 Haumea 系 统 的 运动 来 探究 其 光 心 轨迹 的 
变化 规律 . 仿真 函数 需要 设置 主星 和 卫星 中 心 位 
置 、 大 气 视 宁 度 及 光 通 量 峰值 等 参数 ,以 得 到 二 
维 高 斯 点 源 形状 的 仿真 星象 . 其 中 , 大 气 视 宁 度 使 
用 星象 的 半 高 全 宽 (Full Width at Half-Maximum, 
FWHM) 表 示 . 为 了 控制 单一 变量 , 本 实验 暂 不 考 


二 


Shttps://photutils.readthedocs.io/en/stable/api/photutils.datasets.make_gaussian_sources_ image.html/ 
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虑 添加 泊 松 噪声 4146. 将 主星 置 于 图 像 中 央 ， 从 IM- 
CCE 历 表 中 以 1 ph 为 间隔 获取 从 2022 年 2 月 7 日 到 3 
月 30 日 共 52 d 约 1 个 周期 的 之 卫星 Hi'iaka 相 对 于 主 
星 的 位 置 偏差 ,并 以 紫金 山 天 文 台 姚 安 观 测 站 
0.8 m 和 望远镜 CCD 比 例 尺 0.348(”)-pixel" :换算 出 


卫星 的 像素 坐标 . 根据 文献 [16-17|， 当 星象 
的 FWHM 大 于 2 个 像素 时 , 星象 不 欠 采 样 . 即 在 本 


次 仿真 中 视 宁 度 大 于 0.7” 时 , 星象 不 出 现 欠 采样 情 
Uh. 


图 1 主 卫星 系统 运动 规律 图 . 主星 P 和 卫星 S 同 时 绕 质 心 Cv 做 同一 方向 的 开 普 勒 椭圆 运动 , 光 心 Ci 绕 质 心 Cm 也 做 相同 方向 的 椭圆 运动 . rr 和 rm 
分 别 表示 光 心 和 质心 到 主星 的 角 距 离 , 人 表示 光 心 到 质心 的 角 距 离 , a 表示 主 卫 星之 间 的 角 距 离 . 


Fig.1 Diagram of primary and satellite system. Both the primary P and the satellite S move around the center of mass Cm 


clockwise. The photocenter Cy also moves around Cm clockwise. Symbols rg and rm show the angular distances from the 


photocenter and the center of mass to the primary, respectively. Symbol A shows the angular distance between the photocenter 


and the center of mass. Symbol a shows the separation between the primary and the satellite. 


7£1.0". 1.5’. 2.0". 2.5”. 3.0K Ml rE TF, 
根据 主 卫 星 平均 星 等 差 2.76 mag 设 置 主 卫星 的 光 
通 量 , 模拟 亮 卫星 Hi'iaka 绕 主星 运动 的 过 程 , 并 采 
用 二 维 高 斯 函数 (Gaussian，GS)H6 对 混合 星象 拟 
合 得 到 光 心 位 置 , 结果 如 图 2 (a) 所 示 . 视 宁 度 (用 星 
象 FEWHM 表 示 ) 数 值 越 大 , 光 心 轨迹 形状 越 接 近 于 
WALL, 但 其 大 小 仍 小 于 真实 的 光 心 轨迹 (粉色 ). 并 
且 当 视 宁 度 小 于 2.0/ 时 其 明显 同 中 心 对 称 凹 陷 , 这 
与 2.1 节 所 述 的 光 心 轨迹 规律 四 并 不 相符 . 该 现象 可 
能 与 所 用 的 定 心 算法 有 关 . 当 卫 星 距 离 主星 最 远 
时 , 混合 星象 出 现形 状 不 规则 的 情况 , 且 视 宁 度 数 
值 越 小 越 不 规则 . 因此 使 用 二 维 高 斯 定 心 算法 进行 
定 心 会 受到 一 定 的 影响 . 

为 了 验证 这 种 猜测 , 本 文 也 使 用 了 修正 窍 算法 
(Modified Moment, MM)L9 确 定 光 心 . 修正 矩 算法 
通过 求 一 阶 算 确定 目标 的 质心 , 目标 的 光度 质心 就 
是 光 心 , 两 者 的 定义 相符 . 但 其 对 于 星象 靠近 拟 合 


区 域 边 界 剖 来 的 位 置 误差 最 为 敏感 "9, 所 以 采取 
PR UR AES TE FIERA Ta SK, 即将 第 一 次 的 初步 定 心 结 
果 作 为 第 二 次 的 输入 , 来 保证 星象 位 于 拟 合 区 域 中 
央 . 采用 修正 秆 算法 定 心 并 重复 上 述 过 程 , 光 心 轨 
迹 仿真 结果 如 图 2 (b) 所 示 . 

显然 , 修正 矩 定 心 下 的 光 心 轨迹 在 不 同 视 宁 度 
时 相同 并 且 和 真实 的 光 心 轨迹 (粉色 ) 重 合 . 这 表明 
{ETE FE rE Le 2 RIF AN SSL TER Me, 准确 测量 
TJEME. (IEA RIAR IAA i SUS 
区 域内 像素 到 质心 距离 的 一 次 方 成 正比 , 此 线性 关 
系 使 得 星象 的 FWHM 大 小 变化 时 ( 即 视 宁 度 变化 )， 
其 质心 位 置 不 变 . 

下 面 将 对 这 两 种 定 心 方法 确定 的 光 心 位 置 进 
行 定量 分 析 . 每 种 定 心 方法 都 仿真 了 5 个 不 同 视 宁 
度数 值 下 的 卫星 绕 转 周期 (每 周期 有 1248 个 光 心 位 
置 点 ), 计算 出 光 心 到 主星 的 角 距 离 7L 作 为 纵 坐 标 ， 
主星 和 卫星 之 间 的 角 距 离 a 作 为 模 坐 标 , 如 图 3 所 
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示 . WEIE RB EL (MM) 下 7L 与 之 间 的 关系 在 不 同 视 
宁 度 下 均 为 相同 的 直线 , 且 7i 实 测 值 (MM) 因 与 2.1 
节 (1) 式 计算 得 到 的 真实 值 (real) 相 同 而 重合 . 高 斯 
定 心 (GS) 下 六 与 a 之 间 的 关系 在 不 同 视 宁 度 下 整体 
呈现 不 同 倾 角 的 曲线 , 此 定 心 结果 明显 受到 了 大 气 
视 宁 度 变 化 的 影响 . 当 卫 星 距 离 主星 越 远 及 视 宁 度 


satellite point 

photocenter: FWHM=1.0" 
photocenter: FWHM=1.5 " 
photocenter: FWHM=2.0 " 
photocenter: FWHM=2.5 " 
photocenter: FWHM=3.0" 
real 


© 
Nt 
\ 


\ 


x / pixel 


(a) GS 方法 
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数值 越 小 时 ,xi 实测 值 与 真实 值 相差 越 大 . 在 本 次 
仿真 中 高 斯 定 心 (GS) 下 的 立 实 测 值 与 真实 值 相差 
最 大 为 ~0.284 pixel (0.099), 其 受到 视 宁 度 变化 的 
影响 最 大 为 ~0.214 pixel (0.074”). 综合 定性 和 定量 
分 析 , 在 无 噪声 干扰 时 修正 算 方 法 比 二 维 高 斯 方法 
更 适合 测量 主 卫星 系统 的 光 心 位 置 . 


satellite point 

photocenter: FWHM=1.0" 
photocenter: FWHM=1.5 " 
photocenter: FWHM=2.0" 
photocenter: FWHM=72.5 " 
photocenter: FWHM=3.0" 
real 


x / pixel 


(b) MM 方 法 


图 2 左 、 右 两 幅 子 图 分 别 展示 了 高 斯 方法 和 修正 矩 方 法 定 心 下 的 光 心 轨迹 . 其 中 , 紫色 轨迹 为 仿真 时 卫星 绕 转 主星 一 周 的 轨迹 , 粉色 轨迹 为 真实 的 
光 心 轨迹 . 右 下 角 为 红色 矩形 内 部 的 细节 信息 , 图 (a) 中 从 里 到 外 依次 是 视 宁 度 (FWHM) 数 值 为 1.0”、1.5”“、2.0”“、2.5”“、3.0” 时 的 光 心 轨迹 ， 
图 (b) 中 不 同 视 宁 度 下 的 光 心 轨迹 因 相 同 而 互相 重 薛 . 


Fig.2 The left and right panels show the photocenter tracks with GS and MM centering algorithms, respectively. In each panel, 


the purple track shows the orbit of the satellite around the primary during the simulation, and the pink track shows the real 


photocenter track. The details of the interior of the red rectangle are shown in the lower right inset. The photocenter tracks in 
panel (a) are shown from the inner to exterior with the seeings (FWHMs) of 1.0”, 1.5”, 2.0”, 2.5”, 3.0”, and the photocenter 


tracks in panel (b) are the same at different seeings and overlaps. 


2.3 ”本 文 定 心 算法 概述 

本 文 同样 仿真 了 有 噪声 的 情况 . 假定 主星 与 卫 
星 只 在 图 像 z 轴 方向 相 隅 4 个 像素 (相当 于 IMCCE 
历 表 中 主 卫 星 的 最 远 角 距离 1.45”). 在 该 情况 下 , 两 
种 定 心 算法 下 的 光 心 位 置 差 异 最 大 . 本 文 参考 了 实 
际 观测 图 像 的 背景 , 对 每 幅 仿真 图 像 添加 了 背景 噪 
声 以 及 星象 噪声 , 均 为 随机 泊 松 噪声 . 增加 噪声 后 
目标 星象 的 信 噪 比 为 ~35, 与 姚 安 观测 站 0.8 m 望 远 


镜 观 测 Haumea 系 统 的 信 噪 比 相 近 . 视 宁 度 在 1”-3” 
间 以 0.01” 为 间隔 顺序 取 值 . 在 每 一 视 宁 度 下 仿真 
无 噪声 的 主 卫 星系 统 并 用 高 斯 和 修正 算 方 法 分 别 
定 心 , 接着 添加 20 次 的 随机 泊 松 噪声 并 再 次 使 用 上 
述 两 种 方法 定 心 , 计算 星象 z 轴 坐标 的 均值 和 标准 
Fe. 将 有 无 噪声 时 z 轴 坐标 的 偏差 作为 纵 坐 标 , T 
度 作 为 横 坐 标 , 探究 噪声 对 高 斯 和 修正 矩 方法 定 心 
结果 准确 性 的 影响 , 结果 如 图 4 所 示 . 高 斯 定 心 受 到 
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噪声 的 平均 影响 为 ~1 mas, T4 TE 4E xe Cy BS Bl] 
声 的 平均 影 啊 高 达 ~38 mas. 


m FWHM: 1.0 "(GS) 
FWHM: 1.5 "(GS) 
FWHM: 2.0 "(GS) 
FWHM: 2.5 "(GS) 

@ FWHM:3.0"(GS) 


ri / pixel 


2.5 3.0 3.5 


a / pixel 


图 3 光 心 到 主星 的 角 距 离 rL 与 主 卫 星 角 距离 a 之 间 的 关系 


Fig.3 The relationship between the angular distance ry, 
from the photocenter to the primary and the angular 


distance a from the satellite to the primary 


offset in x /" 


0.001" + 0.010" 


1.5 2.0 2.5 3.0 


seeing /" 


图 4 BWA EED AE OS BRE BY Fe PE 


Fig.4 The influence of noise on the centering results of GS 


and MM centering algorithms 


因此 , 在 实际 拍摄 的 CCD 图 像 中 , 直接 对 单一 
的 混合 星象 (不 可 分 辨 主 卫星 系统 ) 改 用 修正 矩 定 心 
并 不 能 达到 预期 . 基于 高 斯 定 心 受 噪声 影响 小 而 修 
正德 定 心 在 无 噪声 时 更 准确 的 特点 , 本 文 提 出 新 的 
定 心 算法 (GS_new) 作 交心 校正 , 以 减 小 视 宁 度 变化 
对 光 心 位 置 测 量 市 来 的 影响 . 假设 有 噪声 时 星象 的 
Ty rae ty NGS), 无 噪声 时 星象 的 修正 秆 定 心 和 
高 斯 定 心 分 别 为 MM(YD 和 GSWY, 计算 两 者 间 的 坐标 
Fe, 将 此 坐标 差 补 偿 到 GS42 E, 得 到 校正 后 的 新 坐 
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ASC, 如 (3) 式 所 示 : 
fase) = Gs® + (MM® — Gs), a 


Gs) = GSP + (MMP — Gs). 


3 ”实测 资料 光 心 校正 结果 
3.1 ”观测 资料 及 处 理 过 程 

我 们 从 2022 年 2 月 7 日 至 5 月 25 日 用 3 台 不 同 的 
地 基 光 学 望远镜 观测 Haumea 系 统 , 历时 29 晚 (其 
中 4 上 晚 由 2 台 在 不 同 站 点 的 望远镜 同时 进行 观测 ) 共 
招 摄 463 幅 CCD 图 像 . 所 使 用 仪器 的 相关 参数 如 
表 1. 观测 过 程 中 均 采 用 Cousins-I 滤 光 厂 , 以 减少 闫 
色 较 差 折射 的 影响 al, 但 拍摄 的 CCD 图 像 会 出 现 
FERAT, 在 对 观测 图 像 进 行 了 本 底 、 平 场 、 暗 
场 校正 后 , 按照 文献 18] 的 方法 去 除 干涉 条 纹 . 接着 
通过 半 届 览 法 测算 每 一 幅 处 理 后 CCD 图 像 的 大 气 
视 宁 度 范围 2 , 经 统计 其 分 布 在 1.0”-3.0” 之 间 . 

经 过 预 处 理 , 本 文 对 观测 资料 进行 初步 归 算 . 
首先 采用 二 维 高 斯 定 心 方 法 确定 各 CCD 图 像 所 有 
星象 中 心 的 像素 坐标 , 并 将 恒星 的 像素 坐标 与 Gaia 
Data Release 3 (Gaia DR3)P! Æ X PL Ht. 其 次 计 
算 匹 配 恒星 在 拍摄 时 刻 的 天 体 测 量 位 置 , 采用 加 
权 30 常 数 模型 2 23 求 解 底片 常数 . 根据 底片 常数 ， 
Haumea 系 统 光 心 的 像 了 素 坐标 (z,g 妇 转换 为 标准 坐 
ERCE N) (和 7 分 别 为 亦 经 和 赤 续 在 切 点 处 的 投影 
平面 坐标 骂 )， 再 通过 反 投 影 公 式 约 计算 出 光 心 
的 赤道 坐标 (a, 6) (以 春分 点 为 原点 的 亦 道 坐标 系 ， 
oa 表示 赤 经 ,6 表示 赤 纬 ), 将 其 与 IMCCE 历 表 中 的 
理论 位 置 ( 即 Haumea 系 统 质 心 位 置 ) 比 较 . 最 终 得 
到 此 系统 光 心 的 位 置 残 差 O — C (The Observed 
minus Computed), 即 观测 系统 光 心 与 理论 系统 质 
心 的 偏差 . 在 同一 上 晚 对 目标 的 多 次 观测 中 ,本文 
计算 了 杰 经 和 赤 纬 两 个 方 同 位 置 残 差 的 均值 和 
by HE Ze. 对 于 每 幅 图 像 ，Haumea 的 位 置 测量 误差 
在 20-50 mas 之 间 , 此 误差 具体 取决 于 观测 条 件 . 
3.2 ”主星 瞬时 光 变 对 位 置 的 最 大 影响 

据 文 献 [12], 亮 卫 星 Hi'iaka 的 光度 只 占 主星 光 
ERA. 其 光 变 的 振幅 为 19% + 1%), 根据 亮 卫 
星 Hiiaka 的 光度 占 比 和 光 变 幅度 , 可 估计 出 其 对 系 
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统 光 变 的 影响 小 于 1%. 因此 可 以 认为 系统 光 变 主 量 比 并 进行 MM 定 心 , 与 按照 平均 星 等 其 数值 设 定 
要 是 由 主星 光 变 引起 的 . EP SEO) 对 本 批 资料 进行 主 卫 星 通 量 比 并 GS 定 心 的 结果 相 减 , 得 到 了 补偿 
了 光度 测量 , 得 到 了 Haumea 系 统 的 光 变 曲线 (由 主 值 . 因此 , 补偿 值 同 时 包含 了 视 宁 度 变 化 和 主星 光 
星 引 起 的 光 变 ). 因此 , 根据 光 变 曲线 、 观 测 资 料 的 度 变 化 对 交心 位 置 的 影响 . 经 实验 测量 , 主星 Hau- 
时 间 相 位 和 主 卫星 的 平均 星 等 差 2.76 magl, 可 以 mea 光 变 对 交心 的 位 置 测量 影响 最 大 不 超过 9 mas, 
计算 出 观测 时 刻 主 卫星 的 星 等 差 . 仿真 Haumea 系 与 李 丹 等 2 通过 理论 公式 推导 出 的 数值 一 致 

统 星 象 时 , 按照 计算 的 星 等 差 数 值 来 设 定 主 卫星 通 


表 1 望远镜 及 CCD 的 详细 说 明 
Table 1 Specifications of the telescope and its CCD 


Parameters 2.4 m (Lijiang) 1 m (Kunming) 0.8 m (Yaoan) 
Focal Length 1920 cm 1330 cm 800 cm 
Diameter of primary mirror 240 cm 100 cm 80 cm 
Size of CCD array 1900 x 1900 4096 x 4112 2048 x 2048 
Size of pixel 13.5 um x 13.5 um 15 um x 15 um 13.5 um x 13.5 um 
Scale factor 0.283 (")-pixel™* 0.233 ()- pixel™’ 0.348 (’) - pixel! 
Field of view 9.0’ x 9.0 16.0’ x 16.0 11.8’ x 11.8’ 
Exposure time 100 s 500 s 300-600 s 
Number of images 60 82 321 
Night (s) observed 1 7 25 
3.3 ”结果 到 了 在 赤 经 和 赤 纬 方 同 的 位 置 残 产 结果 , 分 别 如 


将 2.3 节 的 新 定 心 算法 应 用 于 实际 观测 ”图 5 和 图 6 所 示 . 图 中 , 空心 矩形 是 每 晚 观测 资料 初 
的 CCD 图 像 中 , 定 心 的 准确 性 与 仿真 的 真实 性 。” 步 归 算 得 到 的 平均 9 - C, 实心 圆 形 是 经 过 光 心 
密切 相关 . 针对 每 一 幅 CCD 图 像 , 尽 可 能 模拟 实 。 ”校正 后 的 归 算 结果 . 根据 Peng 等 图 的 方法 对 这 两 组 
际 观测 中 的 Haumea 系 统 成 像 . 具体 来 说 , 第 一 步 。 ”位 置 残 差 进行 拟 合 得 到 的 三 参数 结果 如 表 2 所 示 ， 
从 IMCCE 历 表 中 获取 拍摄 时 刻 卫 星相 对 于 主星 第 2 列 为 (2) 式 中 用 于 计算 光 心 到 质心 角 距 离 和 的 
EREMAN AERA), 青 通 过 饮 步 。 比例 系数 , 第 3 列 4 和 第 4 列 D 分 别 为 赤 经 和 赤 纬 方 
归 算 得 到 的 底片 常数 模型 计算 出 CCD 的 取 癌 ( 即 向 上 观测 时 刻 的 真实 质心 相对 于 IMCCE 历 表 理 
图 像 4 轴 与 4 方向 的 夹 角 ) 和 比例 尽 ， 以 确定 主星 。 论 质心 的 偏 移 量 . 将 IMCCE 历 表 中 对 应 时 刻 的 主 
与 卫星 的 中 心 位 置 ; 第 二 步 从 IMCCE 历 表 中 获取 。。 卫星 间距 a 及 表 2 中 每 一 行 的 3 个 参数 代入 文献 四 中 


拍摄 时 刻 Haumea 系 统 在 未 经 和 赤 纬 方 同 的 角 速 的 (3) 式 ， 可 方便 得 到 校正 前 后 位 置 残 差 的 拟 合 线 : 
移动 距离 ; 第 三 步 计 算 拍 摄 图 像 时 的 视 宁 度 ; 第 法 的 拟 合 结果 

四 步 结合 Haumea 系 统 光 变 曲线 26 得 到 主星 光度 — 

变化 下 的 主 卫 星光 通 量 比 ; 第 五 步 联合 上 述 参 A 讨论 

数 , 模拟 CCD 在 曝光 时 间 内 对 Haumea 系 统 的 采样 

过 程 ; 最 后 采用 2.3 节 的 新 定 心 算 法 确定 光 心 , 得 4l ELAR 

到 新 的 像素 坐标 . 以 3.1 节 步骤 重新 进行 归 算 , 得 本 文 使 用 新 定 心 算法 对 原来 的 GS 定 心 结 
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进行 了 校正 以 减少 视 宁 度 和 主星 光度 变化 对 
系统 光 心 准确 位 置 测量 的 影响 . 与 原来 的 GS 定 
心算 法 相 比 , 本 文 提 出 的 定 心 算法 能 更 准确 
地 测量 出 Haumea 系 统 的 光 心 , 具体 体现 在 两 方 
面 . 第 一 , 模型 的 拟 合 优 度量 化 指标 一 归 一 化 x? 
( 表 2 第 5 列 Reduced_x*) 在 本 文 提 出 的 定 心 算法 下 
为 3.28, 要 小 于 使 用 GS 定 心 算法 时 的 结果 3.90. 本 
文 的 方法 在 一 定 程度 上 减 小 了 归 一 化 x* 的 数值 , 拟 
合 效 果 有 一 定 的 提升 . 因为 光 心 位 置 测量 同时 也 受 
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其 他 因素 的 影响 , 例如: Haumea 光 环 、 较 差 颜 色 折 
WAM Kamm, 所 以 归 一 化 x? 值 没有 达到 理 
想 情 况 下 的 数值 1. 本 文 作 者 将 在 未 来 考虑 更 多 因 
素 以 进一步 提高 模型 与 矮 行 星 Haumea 系 统 的 契合 
度 ; 第 二 , 使 用 本 文 提 出 的 定 心 算法 后 , 拟 合 输出 的 
平均 比例 系数 从 0.0475 增 加 a 到 0.0681, 与 原 设 定 参 
数 ( 即 通 量 比 12.71, 质量 比 223.8) 代 入 (2) 式 计算 出 


的 理论 值 0.0684 更 为 接近 . 该 数值 也 表明 了 本 文 定 
心 结 果 准 确 性 的 提升 . 


0.08 
~ 0.04 
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JD-2459000 
图 5 Haumea 在 赤 经 方向 上 的 残 差分 布 图. 空心 矩形 表示 使 用 高 斯 定 心 算 法 时 的 平均 位 置 残 差 并 用 点 线 拟 合 , 实心 圆 形 表示 使 用 新 定 心 算 法 时 的 平 
均 位 置 残 差 并 用 实 线 拟 合 . 


Fig.5 The positional residual distributions of Haumea in right ascension. The hollow rectangles show the average position 
residuals fitted with a dotted line when using the Gaussian centering algorithm, and the solid circles show the position residuals 


fitted with a solid line when using the new centering algorithm. 
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图 6 Haumea 在 赤 纬 方向 上 的 残 差 分 布 图 . 空心 矩形 表示 使 用 高 斯 定 心 算 法 时 的 平均 位 置 残 差 并 用 虚线 拟 合 , 实心 圆 形 表示 使 用 新 定 心 算法 时 的 平 
均 位 置 残 差 并 用 实 线 拟 合 . 


Fig.6 The positional residual distributions of Haumea in declination. The hollow rectangles show the average position residuals 
fitted with a dotted line when using the Gaussian centering algorithm, and the solid circles show the position residuals fitted 


with a solid line when using the new centering algorithm. 
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表 2 三 参数 模型 的 拟 合 结果 
‘Table 2 The fitted result of the three parameters model 
Method k AJ” D/” Reduced_y? 
GS 0.0475 + 0.0012 —0.0271 + 0.0009 —0.0083 + 0.0008 3.90 
GS_new 0.0681 + 0.0011 —0.0288 + 0.0009 一 0.0093 + 0.0008 3.28 


卫星 Hi'iaka 在 JPE 历 表 与 IMCCE 历 表 的 
位 置 对 比 
在 第 3.3 节 拟 合 位 置 残 差 的 过 程 中 , 需要 从 历 
表 获 取 卫 星 Hi"iaka 相 对 于 主星 在 赤 经 和 炙 纬 方 问 
的 角 距 离 参 数 , 本 文 发 现 同 一 时 刻下 赤 经 方向 上 
的 该 参数 在 JPL#110 和 IMCCE 历 表 中 存在 较 大 的 
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图 7 


差异 , 如 图 7 所 示 . 虚线 和 实 线 分 别 是 基于 JPL#110 
与 IMCCE 历 表 得 到 的 拟 合 线 , 两 条 线 几 乎 呈现 完 
全 相反 的 波动 . 显然 , 经 过 交心 校正 后 的 观测 数据 
(表示 炙 经 方 同 的 位 置 残 差 ) 波 动 周 期 规律 与 实 线 
(数据 源 于 IMCCE 历 表 ) 理 论 周期 规律 基本 一 致 . 
此 , JPLX#A110 历 表 中 卫星 Hi'iaka 相 对 于 主星 在 赤 经 
方 回 的 角 距 离 可 能 不 准确 . 


@ GS new (This work) 
fitted_IMCCE 
fitted JPL 


680 700 720 740 
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卫星 Hi"iaka 在 JPL 历 表 与 IMCCE 历 表 中 的 位 置 比 较 . 实心 圆 形 表示 使 用 新 定 心 算法 时 在 赤 经 方向 的 平均 位 置 残 差 . 点 线 和 实 线 分 别 是 基 


于 JPL 与 1MCCE 历 表 得 到 的 拟 合 线 . 


Fig.7 The position comparison of the satellite Hi’iaka in the JPL ephemeris and IMCCE ephemeris. The solid circles show the 


average position residuals in right ascension when using the new centering algorithm. The dotted and solid lines are the fitted 


lines obtained based on JPL and IMCCE ephemeris, respectively. 


总 结 与 展望 

本 文 使 用 仿真 法 模拟 主 卫 星系 统 的 运动 ,以 
柯 伊 伯 市 中 典型 的 不 可 分 辨 Haumea 主 卫星 系统 
为 例 , 探究 了 其 光 心 的 准确 位 置 测量 方法 . 根据 无 
噪声 下 修正 怎 比 二 维 高 斯 方法 能 够 更 准确 地 测量 
系统 光 心 及 二 维 高 斯 定 心 结果 受 噪 声 影响 小 的 特 
点 , 提出 了 一 种 新 的 定 心 方法 并 应 用 于 29 晚 目标 为 
矮 行 星 Haumea 系 统 (Haumea-Hi'iaka) 的 观测 图 像 . 
应 用 结果 表明 , 该 方法 减少 了 视 宁 度 和 主星 光度 变 


5 


化 对 主 卫 星系 统 光 心 位 置 准确 测量 的 影响 , 有 利于 
得 到 其 准确 的 轨道 参数 . 该 方法 也 可 扩展 应 用 于 类 
似 的 主 卫星 系统 或 双星 系统 . 
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Research on the Centering Algorithm of Dwarf Planet Haumea 
System 


HAO Jia-nan'? PENG Qing-yu'? GUO Bi-feng’? CAO Jiang-long’? HONG Yao-yaot? 


(1 Department of Computer Science, Jinan University, Guangzhou 510682) 
(2 Sino-French Joint Laboratory for Astrometry, Dynamics and Space Science, Jinan University, Guangzhou 510682) 


ABSTRACT Affected by the Earth’s atmosphere, the image of primary and satellite system observed by 
ground-based telescopes often appears unresolved. ‘Therefore, there may be some deviation between the 
system photocenter and the center of mass. The accurate measurement of the center of mass is helpful 
to improve the orbital parameters and to reveal the formation and evolution of the solar system. Taking 
the dwarf planet Haumea and its brighter satellite Hi’iaka as an example, we simulate its motion of the 
photocenter around the center of mass and explore whether the track of photocenter is the same with 
different seeings. The simulation results show that when using the two-dimensional Gaussian centering 
algorithm, the photocenter track is changing with the change of seeing, however, this phenomenon will 
not occur when using the modified moment centering algorithm. Based on the simulation results and 
the influence of noise, a new centering algorithm is proposed in this paper, which can effectively reduce 
the influence of changing seeing on the accurate photocenter position measurement. In addition, we also 
consider the influence of primary’s changing luminosity on the photocenter position measurement. ‘The 
images of dwarf planet Haumea were observed over 29 nights from February 7, 2022 to May 25, 2022 by 
three telescopes located in Yunnan Province. Our results show that the fitted error of the new centering 
algorithm’s results is smaller than the two-dimensional Gaussian algorithm’s results. In addition, we find 
that theoretical position of the bright satellite from Jet Propulsion Laboratory (JPL) and Institut de 
Mécanique Céleste et de Calcul des Ephémérides (IMCCE) ephemeris has a large deviation. 


Key words astrometry, techniques: image processing, methods: data analysis, Kuiper belt objects: indi- 
vidual: Haumea 
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